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Abstract. We apply a parameterization-independent approach to fitting the dark 
energy equation-of-state (EOS). Utilizing the latest type la supernova data, combined 
with results from the cosmic microwave background and baryon acoustic oscillations, 
we find that the dark energy is consistent with a cosmological constant. We establish 
independent estimates of the evolution of the dark energy EOS by diagonalizing the 
covariance matrix. We find three independent constraints, which taken together imply 
that the equation of state is more negative than -0.2 at the 68% confidence level in 
the redshift range < z < 1.8, independent of the flat universe assumption. Our 
estimates argue against previous claims of dark energy "metamorphosis," where the 
EOS was found to be strongly varying at low redshifts. Our results arc inconsistent 
with extreme models of dynamical dark energy, both in the form of "freezing" models 
where the dark energy EOS begins with a value greater than -0.2 at z > 1.2 and 
rolls to a value of -1 today, and "thawing" models where the EOS is near -1 at 
high redshifts, but rapidly evolves to values greater than -0.85 at z < 0.2. Finally, 
we propose a parameterization-independent figure-of- merit, to help assess the ability 
of future surveys to constrain dark energy. While previous figures-of-merit presume 
specific dark energy parameterizations, we suggest a binning approach to evaluate dark 
energy constraints with a minimum number of assumptions. 
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1. Introduction 

Distance estimates to Type la supernovae (SNe) are currently a preferred probe of the 
expansion history of the Universe [I] , and have led to the discovery that the expansion is 
accelerating [21 El HI IS]- It is now believed that a mysterious dark energy component, with 
an energy density ~70% of the total energy density of the universe, is responsible for the 
accelerated expansion [6j [7J . While the presence of acceleration is now well established 
by various cosmological probes, the underlying physics remains a complete mystery [8]. 
As the precise nature of the dark energy has profound implications, understanding its 
properties is one of the biggest challenges of modern physics. 

With the advent of large surveys for Type la supernovae, such as the Supernova 
Legacy Survey (SNLsfl [9] and Essence^] [I], among others, it is now possible to consider 
detailed studies of the expansion history of the Universe, and shed light on the underlying 
physics responsible for the acceleration. Although the dark energy may be complex, thus 
far it is generally described by a cosmological constant, or through a simple dynamical 
component such as a single scalar field rolling down a potential [10J. The observational 
data is then used to constrain these simple models, generally in the from of determining 
a dark energy equation-of-state (EOS) describing the ratio of its pressure to its density 

or by measuring dynamical parameters such as the cosmic jerk [12]. Using the 
EOS as the primary variable, several studies have considered how current and future 
data might be used to make statements on the physics responsible for dark energy |13j . 
including attempts to establish the shape of the scalar field potential [HI [15] . 

When model fitting data it is generally assumed that the dark energy EOS as 
a function of redshift, w(z), follows a certain predetermined evolutionary history. 
Common parameterizations include a linear variations with redshift, w(z) = wq + w z z 
[IB] , an evolution that asymptotes to a constant w at high redshift, w(a) = Wo+w a (l— a) 
with a as the scale factor [IT], or an evolution with an EOS of the form w(z) = 
Wq — aln(l + z) [18]. Unfortunately, fitting data to an assumed functional form leads 
to possible biases in statements one makes about the dark energy and its evolution, 
especially if the true behavior of the dark energy EOS differs significantly from the 
assumed functional form [15]. Moreover, statements related to the dark energy EOS 
are often made under the assumption of a spatially flat universe, while there still exists 
percent-level uncertainties on the curvature. 

Instead of using a parameterized form for w(z), one can also utilize a variant of the 
principal component analysis advocated in Ref. [19] to establish the EOS with redshift. 
This was first applied in Ref. [2Q] to a set of supernova data from Ref. p]. Recently, Riess 
et al. [3] analyzed a new set of z > 1 SNe from the Hubble Space Telescope combined 
with low redshift SNe, and by making use of the same technique of decorrelating the Wi 
binned estimates, established that the EOS is not strongly evolving. Since the analysis 
of Riess et al. [3], the supernovae sample size has increased by at least by a factor of 2 

| http:/ /www. cfht.hawaii.edu/SNLS/ 
§ http://www.ctio.noao.edu/wsne/ 
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through the Essence survey [I]. The SNe light curves from several independent datasets 
have been analyzed with a common method to extract distance moduli in Ref. [5], and 
we use this publicly available data[jj] to extract the EOS. 

We combine the different distance measurements to extract independent estimates 
of the EOS when binned in several redshift bins in the redshift range < z < 1.8. 
We make use of these binned and uncorrelated estimates of the EOS to address a 
simple question: is the dark energy consistent with a cosmological constant? For a 
cosmological constant w(z) = — 1 exactly, while dynamical dark energy models lead to 
an EOS that either evolves from a large value to -1 today (models which are categorized 
as "freezing" in Ref. [21]) or from -1 at high redshift to large values today ("thawing" 
models of Ref. |21j). In our analysis we also allow departures from a flat universe, but 
allow curvature to be constrained based on complementary information from the cosmic 
microwave background ( WMAP; [7]) and baryon acoustic oscillation distance scales [26J. 
As discussed in Ref. [20], our analysis is facilitated by the fact that our measurements 
of the EOS are completely uncorrelated. Although we focus on the EOS, we note that 
one can also extract uncorrelated estimates of other parameters related to the expansion 
history and dark energy, such as the dark energy density [22]. However, unlike the case 
with the EOS, such estimates cannot be used to directly address the simple question of 
whether the dark energy is a cosmological constant. 

We make use of our binned estimates to comment on several new developments 
related to dark energy studies. First, in addition to addressing the extent to which 
w(z) = —1 in present data, we also comment on the extent to which dynamical dark 
energy models, such as "freezing" and "thawing" models [21], may be distinguished 
or ruled out with current cdistance data. We also discuss the role uncorrelated, 
independent EOS measurements can play in furthering our understanding of the dark 
energy. Recently, the Dark Energy Task Force (DETF) [23] has suggested a figure-of- 
merit to compare the abilities of different experiments to extract information related to 
dark energy. This is done in terms of the inverse area of the error ellipse of the equation 
of state and its evolution with redshift, utilizing either the Linder parameterization 
[TT] or with two parameters of the form w(a) = w p + (a p — a)w a [13]. The discussion 
is then restricted to a two parameter description of the dark energy equation of state, 
assuming a very specific evolutionary behavior. The ability to extract information about 
dark energy from current and future experiments thus becomes a model dependent 
statement. Using uncorrelated, independent equation of state estimates, we propose 
a model independent figure-of- merit. Our approach involves the inverse of the sum of 
inverse variances of uncorrelated w(zi) bins as an way to capture all of the available 
information related to dark energy in a given data set or experiment. 

The paper is organized as follows: In the next Section we review techniques for 
reconstructing the EOS, we following the methods of Ref. [20] and Ref. [3]. In Section 3 
we present our results, addressing whether dark energy is a cosmological constant or has 

|| http: / /www. ctio.noao.edu / wprojcct /wrcsults / 
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Figure 1. Hubble diagram for type la supernova data used for the present analysis. 
The dataset includes a total of 192 SNe from the recent analysis of light curves in 
Ref. [5]. For comparison, we also separately analyze a subset of 104 SNe within this 
sample that were also used in Riess et al. [3J. The distance moduli used here for this 
subset, however, are different from this original study due to the reanalysis of light 
curves in Davis et al. [5]. The orange line is for our best fit. 



dynamical behavior that leads to an evolution in the EOS. We show that existing data 
rule out extreme forms of dynamical behavior at the 68% confidence level, including 
models that either start with w(z) > —0.4 at z > 1.0 or models that asymptote to 
values of w(z) > —0.85 at z < 0.2. In Section 3.1 we outline a new figure of merit to 
assess the dark energy information content of future experiments, based on uncorrelated 
binned estimates of the EOS. We conclude with a summary of results in Section 4. 

2. Methodology 

A simple way to model the dark component of the universe credited for the accelerating 
expansion is through a modification of the standard cosmological model. We utilize 
the Friedmann equations, and specify the dark energy density and its equation-of-state 
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(EOS). We assume a piecewise constant EOS, with value Wi in each i th redshift bin 
(defined by an upper boundary at zA. We fit the observational data to the luminosity 
distance as a function of redshift. The expression for luminosity distance, cIl(z), depends 
on whether the universe is flat, positively, or negatively curved (i.e., the sign of fife), 
and is given by 
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and F(z) depends on the binning of w(z). For the n redshift bin F(z) has the form 



F(z n > z > z n -i) 
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(3) 



We define the zeroth bin as zq = 0, so the product is unity for redshift z in the first bin. 
For our primary analysis we set z\ = 0.2, Z2 = 0.5, z% = 1.8, and z^ extends beyond 
the surface of last scattering at zqmb = 1089. We assume w(z > 1.8) = —1 and allow 
variation within the remaining three redshift bins. Selecting the cutoff point for z% is 
fairly arbitrary; we found that pushing it back as far as z% = 2.5 does not substantially 
alter the outcome of our analysis. 

In addition to SNe, we also make use of four primary constraints from the literature 
following the analysis of Riess et al.[3], modifying these to account for variations in 
curvature since the analysis in Ref. [3] assumed an a priori flat cosmological model. 
These constraints are: 

• The product of the Hubble parameter h = H /100 and the present local mass 
density fl m from SDSS large scale structure measurements [38], given by fl m h = 
0.213±0.023. In cases where we allow curvature to vary, we either take a flat, broad 
prior in curvature or, to highlight the result under the assumption of a measured 
value for the curvature, we take Qk = —0.014 ±0.017 as derived by WMAP analysis 
by combining WMAP and the Hubble constant [7]. 

• The SDSS luminous red galaxy baryon acoustic oscillation (BAO) distance estimate 
to redshift zbao = 0.35. Here the constraint is on the overall parameter A = 

(Il(z)/(1 + z) is the angular diameter 



CZBAO 



r l^BAO, 



CZBAO 



H E(zbao) 



1/3 



, where r(z) 

distance. The angular correlation function of red galaxies in the SDSS spectroscopy 
survey leads to A = 0.469(^g)" a35 ± 0.017 [26]. Following Riess et al. [3], we use 
the WMAP estimate for the scalar tilt with n = 0.95 [7]. 

The distance to last scattering, at ;?cmb = 1089, written in the dimensionless form 
Rcmb = v /nA J H ° r(z CMB ) where r(z CM B) is the angular diameter distance to the 
CMB last scattering surface. We use the dark energy and curvature independent 
estimate with R = 1.70 ± 0.03 
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The distance ratio between z = 0.35 and last scattering z = 1089 as measured by 
the SDSS BAO analysis [26]: 

1/3 

(4) 



0.35 — " 

r{Z C MB) 

with the value of R . 35 = 0.0979 ± 0.0036. 

When estimating parameters we make use of the y 2 statistic for a particular model 
with parameter set 6 (wi, H , Q m , and in some cases f^): 

JV / theory _ data \ 
n=l \ a « ^ °mt / 

where N is the total number of supernovae in the sample. While our total sample 
includes 192 supernovae, we also extract a subset of 104 supernovae that was analyzed 
previously [3]. While the subsample is for comparison with previous results, our distance 
estimates differ from the original analysis due to a reanalysis of light curves by Davis 
et al. [5] using a common light curve fitting method. When estimating \ 2 we set an 
intrinsic dispersion of order o~- mt ~ 0.1, such that the minimum \ 2 value for the best fit 
model comes to be about one. Using x 2 we calculate the probability P(9) oc exp(— x ^ ) 
and derive constraints by generating Markov chains through a Monte-Carlo algorithm. 
The algorithm generates a set of models whose members appear in the set (or chain) 
a number of times proportional to their likelihood of being a good fit to the observed 
data, after marginalizing over other priors. The likelihood probability functions for 
each independent parameter are generated by simply taking a histogram over the chain. 
We marginalize over H assuming a broad uniform prior over the range [30, 85] km s^ 1 
Mpc" 1 . We also marginalize over Vt m assuming the quoted prior above for Vt m h using 
SDSS large scale structure measurements. 

In the case of dark energy parameters, the redshift binned EOS estimates are 
correlated such that their errors depend upon each other. These correlations in the 
redshift binned EOS estimates, Wi, are captured by the covariance matrix, and this 
matrix can be generated by taking the average over the chain such that: 

C = (ww T ) - (w) (w T ) . (6) 

This covariance matrix is not diagonal as the values found for the various EOS estimates 
are correlated. This is expected, as the integration over low redshift bins in Eq. (TO 
obviously affects the model fit in middle and higher redshift bins. The behavior is also 
such that the best constraints are found for the lowest redshift bin, with higher bins 
having progressively weaker constraints. With the addition of the distance scale to the 
last scattering surface from CMB data, the constraint in the higher redshift bins are 
improved, as seen in prior analyses [3]. 

Instead of discussing w(z) in correlated bins, we follow Huterer & Cooray [20] and 
transform the covariance matrix to decorrelate the EOS estimates. This is achieved 
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by changing the basis through an orthogonal matrix rotation that diagonalizes the 
covariance matrix. We start by the definition of the Fisher matrix, 

F = C" 1 = O t AO (7) 

where the matrix A is the diagonalized covariance for transformed bins. The 
uncorrelated parameters are then defined by the rotation performed by the orthogonal 
rotation matrix q = Ow and has the covariance matrix of A -1 . There is freedom of 
choice in what orthogonal matrix is used to perform this transformation, and we use the 
particular choice that was advocated in Ref. [20J and write the weight transformation 
matrix: 

W = O t A^O (8) 

where the rows are then summed such that the weights from each band add up to unity. 
This choice ensures we have mostly positive contributions across all bands, an intuitively 
pleasing result, so for example we can interpret the weighting matrix as an indication of 
how much a measurement in the third bin is influenced by SNe in the first and second 
bins. We apply the transformation W to each link in the Markov chain to generate a 
set of independent, uncorrelated measures of the EOS and its probability distributions 
as determined by the observables. 

In addition to probability distribution functions for each of the uncorrelated EOS 
estimates, to study the redshift evolution of w(z), we also study the differences between 
these uncorrelated estimates. The probability distribution functions of the difference 
between estimate io, and Wj is generated through P(wdm) oc / P Wi (w)P w .(w + w<as)dw. 
While in the uncorrelated case an estimate u>j is not necessarily associated with a single 
redshift bin, due to support from adjacent redshift bins due to the transformation, any 
significant difference between Wi and Wj can be considered to be evidence for an evolution 
in the dark energy EOS. Moreover, if the dark energy is associated with a cosmological 
constant, then these difference estimates should be precisely zero. 

3. Results 

The results presented in this paper are derived from the statistical analysis of a 
combination of recent supernova surveys, including the Supernova Legacy Survey 
(SNLS) [9], the ESSENCE survey [I], and high-z supernovae discovered by the Hubble 
Space Telescope (HST) [Hj. In particular, we use a total of 192 SNe la measurements 
taken from a combination of supernovae analyzed in Ref. [5] using a common light curve 
fitting method (Figure 1). Also, for comparison, results are presented for the 104 SNe 
la that overlap with the "Gold" data set presented in Riess et al. [3], although the 
distance moduli values we use here for the same subsample are slightly different from 
the values published in the original analysis due to variations in the light curve fitting. 
This analysis includes the four external constraints outlined in the previous section. In 
addition to the standard flat cosmological model generally assumed when making fits 
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Figure 2. Uncorrelated estimates of the dark energy equation of state using 
a combined sample of supernova data and constraints from WMAP and BAO 
measurements. In the top panel a flat universe prior is assumed (fl^ = 0), and the 
filled and open symbols show the constraint with the total sample of 192 SNe from 
Ref. [5] and a subset of 104 SNe corresponding to a previous analysis [3], respectively. 
In the bottom panel we show w(z) estimates without the flat universe prior; open and 
filled symbols showing constraints with a broad prior for fi^ arid = —0.014 ± 0.017, 
respectively. 



to dark energy EOS, we also allow for variations in the curvature, both with a prior on 
fifc based on WMAP and Hubble constant measurements, and with no prior. 

In Figure 2 we highlight our results for w(z), in the redshift bins z < 0.2, 
0.2 < z < 0.5 and 1.8 > z > 0.5, with w(z > 1.8) = —1. These binned estimates 
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Figure 3. Probability distribution functions of uncorrclated estimates of the dark 
energy equation of state (top panel) and the difference between estimates of the 
equation of state (bottom panel). Solid and dashed lines arc for the case with 192 
and 104 SNe (see caption of Figure 1), respectively. 



are uncorrelated following the technique of Huterer & Cooray [20] and as outlined in 
Section II. The procedure to decorrelate binned estimates that are predefined over a 
certain range usually results in adding contributions from nearby bins, but these are 
generally smaller than the main contribution from the bin in which the estimate was 
first defined. On average, the first uncorrelated measure is 73% determined by its own 
bin, with a minimal contribution from the third bin. The second measure, on average, 
was 50% determined by the second bin, with a substantial contribution from the first 
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Figure 4. Probability distribution function of cosmic curvature by combining 
supernova data with additional constraints outlined in the paper. We assume a broad, 
flat prior for flk, but the combination of SN data and CMB and BAO distance scales 
results in a tight constraint on curvature. We find f2fc = 0.004 ± 0.019 (192 SNe) 
and fifc = 0.002 ± 0.022 (104 SNe), which is fully consistent with estimates made by 
the WMAP analysis by combining WMAP data with supernova data from the SNLS 
survey. 



bin. The third measure was typically 54% determined by the third bin, with a 6% 
contribution from the first bin. 

These results do not presume a particular evolutionary history, as opposed to model 
fitting to a specific form, e.g. w(z) = w + W\Z [16] or w{z) = w + w a {l — a). 
Our procedure, which fits binned values of w(zi) and then decorrelates them, has 
the advantage that one can extract redshift evolution independent of a model. This 
is particularly effective if the model to be assumed turns out not to be an accurate 
representation of the true underlying EOS. 

As shown in Figure 2, the dark energy EOS as a function of redshift is fully 
consistent with w(z) = —1 at the la confidence level, for both the full sample of 192 SNe, 
and the subset of 104 SNe corresponding to the earlier analysis [3]. As shown in the lower 
panel of Figure 2, this conclusion is unchanged when we drop the assumption related to a 
flat cosmological model, regardless of the assumed prior on Qk- Our external constraints 
related to distance to the last scattering surface from CMB, and BAO distance scale to 
z = 0.35, provide a strong constraint on Qk- We explicitly include an additional flatness 
constraint to allow comparison with earlier work [5j [3] . To highlight the extent to which 
EOS estimates are consistent with w(z) = —1, in Figure 3 we plot the probability 
distribution functions P(wi) both for flat cosmologies and for the case with curvature 
allowed to vary. Except for the third bin, which still remains mostly undefined with 
a very broad probability function, the first two bins are peaked and allow constraints 




Figure 5. Left: A comparison of w{z) estimates and dynamical dark energy models, 
based on the Monte-Carlo modeling approach of Ref. [28], that are inconsistent with 
current estimates. We show both cases of "thawing" and "freezing" models inconsistent 
with current data (see text for details). Right: The shapes of potentials generally 
corresponding to w(z) models shown in the left panel which are inconsistent with our 
estimates of the binned EOS from a combined sample of supernovae and cosmological 
distance scale measurements. We show V(<p) as a function of the scalar field 4> when 
potentials are normalized to the value at z = 3 (Vo = V(z = 3)), which we take to be 
= 0.2. 

at a high confidence level over —2 < w < 0. These show a clear consistency with 
w — — 1, and hence are completely consistent with a cosmological constant. While the 
allowed range for w(z 3 ) is broad, at the 68% confidence level we find that w(z 3 ) < —0.2, 
suggesting that we can rule out a large EOS even at z > 0.5. 

The uncorrelated binned estimates of w(z) derived in Huterer & Cooray [20] using 
an earlier "Gold" sample from Riess et al. [29] showed an equation of state that 
varied significantly between the lowest-redshift bin and the second bin. This difference 
decreased in the most recent "Gold" sample as analyzed by Riess et al. [3]. In the 
current work, utilizing an extended sample of supernovae, we no longer find evidence 
for a variation in the dark energy EOS between the first and the second bin. To show this 
explicitly, we also plot the probability distribution function of the difference between 
binned estimates of Wi(z) in Figure 3. Between the first and the second bin we find 
the difference to be W2 — W\ = 0.06 ± 0.26 at the 68% confidence level for the 192 
SNe sample with a flat model. Current data is thus completely consistent with a 
cosmological constant. Previous estimates of a large and a statistically significant value 
for W2 — wi, with wi < —1 and wi > —0.8, led to suggestions in the literature for a 
physical mechanism called dark energy "metamorphosis" [30] . While earlier conclusions 
were limited to a small set of supernovae data, with the larger sample it is now clear 
that there is little evidence for a sudden transition in the EOS around z ~ 0.2. Future 
data could tighten these constraints, either further narrowing down to a cosmological 
constant, or providing evidence for small variations in the EOS with redshift. 
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Figure 6. A comparison of w(z) values from dynamical dark energy models, based 
on the Monte-Carlo modeling approach of Ref. [28]. Each dot represents a potential 
dynamical dark energy model. The data points are the allowed range from our analysis. 



We show that our conclusions are generally unchanged by assumptions related to 
the curvature. This is because we constrain the EOS using a combination of supernova 
data and existing measurements of the cosmic distance scale out to z = 1089 and 
z = 0.35 with CMB and BAO, respectively. The combination of supernova data 
and these measurements, combined with our prior on the Hubble constant, leads to a 
strong independent constraint on the curvature parameter Q^. We show the probability 
distribution P(f2&) in Figure 4 for both the full sample and the subset of 104 supernovae. 
In both cases is consistent with zero; with the full supernovae sample, we find 
fi fc = 0.004 ± 0.019. This is about la away from the combined WMAP and Hubble 



Dark Energy Evolution 



13 



constant estimate of = —0.014 ± 0.017 or the combined WMAP+SNLS estimate of 
tt k = -0.011 ± 0.012 [7]. On the other hand, the combined WMAP+SDSS estimate 
from the same analysis is Q K = — 0.0053±g;gg|g, which is a shift in the direction of the 
fl k value we find when the combined SNe dataset is analyzed with the WMAP and BAO 
priors. 

To demonstrate how our estimates of w(zi) can be used to understand the redshift 
evolution of the dark energy component, in Figure 5 we show a sample of predictions 
related to dynamical dark energy models that are ruled out with present estimates of 
Wi(z) at the 68% confidence level. These cases generally involve a dark energy EOS that 
starts as w > —0.2 at high redshifts, or an EOS that stars with a value around -1 at 
high redshift but evolve to a value greater than -0.85 when z < 0.2. The former models 
belong to the general category of "freezing" models described in Caldwell & Linder 
|21j . while the latter models are categorized as "thawing" models. We generate these 
models following the numerical technique of Huterer & Peiris [28], writing the Klein- 
Gordon equation for the evolution related to a scalar field as </> + 3H<p + dV/d<f> = 0, 
and then numerically generating a large number of models following the Monte-Carlo 
flow approach as used for numerical models of inflation [3TJ [32} [33], EU [35] . We do not 
reproduce details as the process is similar to the modeling of Ref . [28j [36] . 

As a further application of our estimates of the EOS, and the difference between 
two binned estimates of EOS, in Figure 6 we compare our measured value with values 
expected for a large number of models. Again, we rule out certain extreme models 
where w(z{) varies rapidly between adjacent bins. In dynamical dark energy models, 
models that generally lead to a large variation in w between two adjacent bins also have 
a value significantly different from -1 in one of the bins. Thus most models are currently 
ruled out from the value of w in a single bin, rather than through the difference of 
W2 — Wi or w% — io 2) since the the latter are still largely uncertain. While we can 
use the numerical modeling technique of Huterer & Peiris [28] to make qualitative 
statements about the EOS, and to rule out extreme possibilities for its dynamical 
evolution, given the stochastic nature of model generation we cannot use this sort of 
method to make detailed statements about, for example, the scalar field (quintessence) 
potential responsible for dark energy. Instead, it is necessary to directly reconstruct 
the scalar field potential from supernovae distance data. While there are attempts to 
recover the potential by directly model fitting various parametric forms of the potential 
as a function of the scalar field, such as power-law or polynomial functions of the scalar 
field <j), model independent binned estimates of the potential are preferable |37] . 

3.1. A New Figure of Merit 

As discussed in the previous sections, our binned estimates allow us to study the redshift 
evolution of the dark energy EOS without the need to assume an underlying model. This 
is to be contrasted with the usual approach, in which a parameterized form for w(z) is 
required to fit the data. With an increasing supernova sample size, and improvement 
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in other cosmological observations, we may be able to recover three or more binned 
values at the 10% level or better. In the context of planning dark energy experiments, 
and assessing the constraining power of future data, it may be advantageous to consider 
binned estimates of the dark energy EOS, rather than the error associated with a specific 
and arbitrary two parameter model of the equation of state. The latter is the approach 
adopted to quote the "figure of merit" (FOM) of an experiment, based on the inverse of 
the area of the ellipse of the two parameters describing the EOS with redshift (introduced 
in [23]). In the case of binned estimates, once uncorrelated, we can quote an alternative 
FOM as \J2 i l/a 2 (w Zi )] 1 , which takes into account the combined inverse variance of 
all independent estimates of the EOS. In an upcoming paper we will quantify the exact 
number of w(z) estimates that can be determined with future experiments involving 
supernovae and large-scale structure (weak lensing, baryon acoustic oscillations over a 
wide range in redshift), and we will compare this alternative FOM to the method of 
Ref. [23]. 

4. Summary 

We use a sample of 192 SNe la (and a subset for comparison) to constrain the dark energy 
equation-of-state parameter and its variation as a function of redshift. We use a model 
independent approach, providing uncorrelated measurements across three redshift bins 
below z = 1.8, and find that w(z) is consistent with a cosmological constant (w(z) = —1) 
at the 68% confidence level. At the same confidence level we find that the EOS is greater 
than —0.2 over the redshift range of < z < 1.8. Overall, there is no strong evidence 
against the assumption of a flat = universe, especially when recent supernova 
data are combined with cosmological distance scale measurements from WMAP and 
BAO experiments. We argue against previous claims in the literature of evolving dark 
energy, such as dark energy "metamorphosis", where the EOS changes significantly 
from w > — 1 at z > 0.2 to w < — 1 at z < 0.2. Instead, we find consistency with a 
cosmological constant, encapsulated in the 68% level constraint: w 2 — Wi = 0.06 ±0.26, 
where W\ is the value of the dark energy EOS in the z < 0.2 bin, and W2 is the value 
in the bin 0.2 < z < 0.5. A transition in the EOS can also be ruled out between our 
second and third binned estimates of EOS, although we still find large uncertainties in 
our determination of the EOS at z > 0.5, and we are insensitive to rapid variations at 
z > 1. We compare our EOS estimates to Monte Carlo generated dynamical dark energy 
models associated with a single scalar field potential. Our EOS estimates generally 
allow us to rule out extreme "thawing" and "freezing" models, though a large number 
of potential shapes remain in agreement with current data. 

We also suggest an alternative, parameter independent figure-of-merit, with which 
to evaluate the potential of future missions to constrain properties of the dark energy. 
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